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RESUMEN 

La aplicacion de varios metodos para obtener resultados confiables con el fin 
de aproximarnos de forma certera al comportamiento de una estructura ante una 
solicitacion dinamica como es un sismo, usando acelerogramas de sismos pasados. 

La finalidad de este artfculo es interpretar de mejor manera como se obtuvo el 
analisis del metodo de Newmark y el metodo de Wilson, para aplicarlos en una 
estructura y observar su comportamiento con la conclusion de comparar los 
resultados. 

Palabras Claves: Metodo de Newmark, Metodo de Wilson, Solicitacion Dinamica, 
Acelerograma. 


SUMMARY 

The application of several methods to obtain reliable results in order to come 
closer of accurate form the behavior of a structure before a dynamic request as it is 
an earthquake, using accelerographs of past earthquakes. 

The purpose of this article is to interpret of better way as there was obtained the 
analysis of the method of Newmark and the method of Wilson, to apply them in a 
structure and to observe its behavior with the conclusion of comparing the results. 
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1. INTRODUCCION 


Uno de los metodos mas utilizados para encontrar la respuesta en el tiempo 
de una estructura es el denominado Metodo /? de Newmark, razon por la cual en 
este artfculo se presenta el marco teorico respective, para que sea utilizado en forma 
apropiada habida cuenta que no siempre es estable el metodo. En efecto, para 
cuando se halla la respuesta con varios sismos puede darse el caso que con uno 
de ellos, se obtienen resultados que no son logicos que no se presentaban con los 
otros sismos. (Meruane 2012, Aref 2016, Chavez 2010). 

En este caso, se puede introducir amortiguamiento al metodo para lograr 
suavizar la respuesta en cierto intervalo de tiempo y con ello obtener resultados 
logicos y confiables. Otra opcion cuando se tienen resultados que no son logicos es 
utilizar otro metodo de analisis y en este contexto puede ser que se decida por el 
Metodo 6 de Wilson, que es muy parecido al p de Newmark, por lo que tambien se 
presenta el desarrollo numerico de este metodo. 

Los metodos explicados en el presente artfculo tienen como finalidad arrojar 
resultados en el tiempo en funcion de aceleraciones establecidas inicialmente, por 
medio de interacciones con ecuaciones deducidas de forma matematica, partiendo 
de varias hipotesis y colocando variables para estabilizar el calculo, evitando 
resultados incoherentes con respecto a la realidad; ayudado de herramientas 
computacionales como es el caso de Mat Lab (Angulo et al. 2001) y programas 
matriciales de calculo como CEINCI-LAB (desarrollado en la Universidad de 
Fuerzas Armadas ESPE) con lo que se obtiene resultados de forma rapida y 
eficiente para ser interpretados con criterio. 

2. MANUAL RAPIDO DE CEINCI-LAB 

En Aguiar (2012, 2014) se encuentra el marco teorico de Analisis Matricial 
de Estructuras y Dinamica de Estructuras con la presentation de programas 
desarrollados en Mat Lab que forman parte del sistema de computation CEINCI- 
LAB, con los cuales se facilita notablemente encontrar la respuesta en el tiempo de 
una estructura con base empotrada o con aisladores sfsmicos o con disipadores de 
energfa. 


Para usar CEINCI-LAB, el usuario debe saber como se halla la respuesta 
en el tiempo y con este conocimiento saber la secuencia de uso de los programas 
de la librerfa de CEINCI-LAB, los datos que debe dar a cada uno de ellos y los 
resultados que se esperan. 

En este apartado se presentan unos pocos programas con los cuales se 
puede elaborar un programa para encontrar la respuesta en el tiempo de la 
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estructura con base empotrada que se analizara ante el sismo de El Centro de 1 940. 

[NI,NJ]=gn_portico(GEN)._ El programa usado con el vector de colocacion 
gen que tiene como datos el numero del elemento, nudo inicial y nudo final, 
incremento de elementos y la numeracion. 

[X,Y]=glinea_portico(NUDOS)._ El programa genera las coordenadas de los 
nudos de forma plana ayudado de la matriz nudos ingresada con el numero del 
nudo, posicion con respecto a un eje de referenda. 

[CG ngl]=cg_sismo(nod,np,nr)._ El programa realiza una matriz que indica 
los grados de libertad en cada nodo tomando en cuenta que inicia con coordenadas 
principales para una analisis sismico y requiere como datos el numero de nodos, el 
numero de pisos y el numero de restricciones o nudos con apoyos. 

Dibujo(X,Y,NI,NJ)._ El programa realiza un grafico o esquema del portico piano 
para revisar si los nudos y elementos fueron colocados de forma correcta, de esta 
forma el analisis continue el proceso. 

[L,seno,coseno]=longitud(X,Y,NI,NJ)._ El programa usa las coordenadas de 
los nudos y la ubicacion de los elementos para obtener el valor de la longitud de 
cada uno y su ubicacion dentro del piano por medio de un angulo con referencia a 
la horizontal. 

[VC]=vc(NI,NJ,CG)._ El programa genera un vector de colocacion para cada 
elemento que consta de los grados de libertad que tiene en el nudo final y los grados 
de libertad del nudo inicial usado para el ensamble de la matriz de rigidez. 

[ELEM]=gelem_portico(SECCION)._ Este programa genera elementos de un 
portico piano usando la matriz seccion que tiene como datos el numero del 
elemento, las dimensiones de la seccion transversal y los elementos a generar. 

[K]=krigidez(ngl,ELEM,L,senos,coseno,VC,E)._ El programa realiza una 
serie de calculos usando los programas anteriores que nos da como resultado las 
rigideces de los elementos usando el modulo de elasticidad del material, estas 
rigideces son ensambladas de forma directa por medio del vector de colocacion, 
posterior a este paso se usa esta matriz para condensarla en una matriz de rigidez 
lateral en dependencia del numero de pisos o coordenadas principales. 

[T,fi,OM]=orden_eig(KL,M)._ El programa calcula los valores y vectores 
propios de los periodos, los modos de vibracion y la frecuencia de cada piso en 
funcion de la matriz de rigidez y la matriz de masa, estos resultados se ordenan de 
mayor a menor. 
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[C]=amortiguamiento(M,fi,OM,zeda)._ El programa determina la matriz de 
amortiguamiento de la estructura usando el metodo de Wilson y Penzien para lo 
cual se necesita de la matriz de masas, los modos de vibracion, las frecuencias de 
vibracion y el factor zeda de amortiguamiento de la estructura que depende del 
material. 

• Base de datos del sismo 

El analisis con el metodo de Newmark y Wilson depende principalmente por 
conocer las aceleraciones de un sismo dado, estos datos ingresan en la primera 
parte del calculo para determinar la fuerza junto con la masa y el vector J para 
analisis piano. 

La base de datos con aceleraciones en funcion del tiempo de un sismo que se 
escoge para este analisis, es el Sismo del Centro en 1940 que tiene un 
acelerograma del suelo referente al siguiente grafico con una aceleracion maxima 
aproximada de 300 gals (cm/s 2 ). 



Figura 1 Acelerograma de la base de datos Sismo del Centro de 1940. 


Los programas del sistema de computation CEINCI-LAB se hallan 
publicados en el portal de la Universidad Laica Eloy Alfaro de Manabf, ULEAM, 
dentro de la carpeta de Investigaciones se encuentra la carpeta Centro de 
Investigaciones en Ingenierfa Slsmica y dentro de esta se halla otra denominada 

CEINCI-LAB 
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3. METODO DE NEWMARK 

El Sistema de ecuaciones diferenciales que es usado en los problemas 
dinamicos define la ecuacion: 

Mq + Cq + Kq = MJa(t) (1) 

Donde: M es la matriz de masa del sistema; K es la matriz de rigidez del 
sistema; C es la Matriz de amortiguamiento del sistema; q es el vector aceleracion; 
q es el vector velocidad; q es el vector desplazamiento; J es el vector de incidencia 
del movimiento del suelo en los grados de libertad, que contiene valores unitarios y 
ceros, depende del modelo de analisis y a (t) Aceleracion de movimiento de suelo. 

Ya que las aceleraciones, velocidades y desplazamientos se determinan de 
forma interactuante con un At es decir una variacion del tiempo la ecuacion dinamica 
queda: 


Mq i+ i + Cq i+1 + Kq i+1 = MJ a i+1 


( 2 ) 


Se determina la ecuacion de variacion de desplazamiento 

<7i+i = &q i+1 + <?i (3) 

Las ecuaciones de aceleracion y velocidad en funcion de At son: 


1 A 1 

^ i+1 “ P At 2Aqi+1 P At^ 



U'i 

(4) 

, '« = 04t A ‘"« + ( 1 p) q, + t 

t 1 " 

h) 

Attji 

(5) 


Donde y,p, son parametros del metodo que se indicaran posteriormente. Las 
ecuaciones de rigidez efectiva y fuerza efectiva se obtienen al reemplazar (3), (4) y 
(5) en (2) y concentrarterminos. 


K * A<7 ;+ i — F i+1 


( 6 ) 


Rigidez efectiva del sistema de analisis. 


K 


K + M 


1 Y 

£At 2 + C Jit 


(7) 


Fuerza efectiva del sistema de analisis. 
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Fi 4-1 = 


. At 4 Y\ . . (4 1 \ . 


( 8 ) 


3.1 Estabilizar el metodo p de Newmark 

Para iniciar con la demostracion se genera un grafico de la aceleracion para 
un periodo de tiempo At. 



Figura 2 Supuesto de aceleracion constante metodo de Newmark 


Ya que el supuesto de la aceleracion es constante para este analisis se 
considera la ecuacion de la recta como un promedio entre las aceleraciones de los 
dos puntos: 


Ri + t — 2 (Qi+bt F 9 i) 


(9) 


Se integra en funcion de r la ecuacion (9) y se reemplaza t = At con lo que se 
obtiene. 


At .. 

4i+At Ri ~ k ~2 \Ri+At Q i) 


( 10 ) 


Se integra dos veces en funcion de r la ecuacion (9) para determinar el 
desplazamiento y se reemplaza r = At con lo que se obtiene. 

At 2 (11) 

r?i+At — r?i F At cfa + ^ (.Qi+At <7 u 


La ecuacion (10) se remplaza valores de Vz = y y en la ecuacion (11) se 
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reemplaza valores de % = p, dando como resultado las siguientes ecuaciones. 

<?;+At =q i + At(yq i+At + (l-x) ( ? i ) (12) 

<?;+At = q t + At q t + At 2 ( /? q i+At + - /?) q f ) 1 1 

Para facilidad de calculo se establece de forma inicial que 0 = w At, donde uj 
es la frecuencia, se inicia con el producto At con la ecuacion (13) como se muestra. 

Atq l+At = At<7i + At 2 (yq i+At + (l-y)<7' i ) (14) 

Se considera para este caso el uso de la ecuacion dinamica sin el 
amortiguamiento y se despeja la aceleracion recordando que: 



entonces 


U ) 2 = M -1 * K 


Mcji + Kq t = F t 

q t = —coqi + M~ 1 F l ; en terminos de i + 1 resulta 
Ri+i = + M _1 F i+1 

Si no se considera la carga y se reemplaza el valor de q i+1 en la ecuacion (14) 
y el valor de q t en la ecuacion (13), se obtiene lo siguiente: 

At q l+1 + 9 2 y q i+ i = A tq t - 0 2 (1 - y)<fr (15) 

Ri+ 1(! + P O 2 ) = <h + At q t - d 2 - /?) q t ^ ^ 

Las ecuaciones mencionadas en la parte anterior se las escribe en forma 
matricial, se coloca una letra para identificarla y agilitar el calculo. 


Xi + pe 2 ) 

01 r di +1 1 _ 

i-0 2 (^-p) 

1 r i 

. e 2 y 

lJ lAt q i+1 ) ~ 

\Z ) 

-e 2 (i -r) 

lA tqj 

[B] 

Xt+i 

[C] 

*t 


X t+1 = [B] _1 [C] X; 


El valor de [£] -1 [C] se le conoce como operador de integration y sera llamado 
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A, se requiere como condition que p(A ) < 1 , entonces el determinante de la matriz 
C menos la matriz B multiplicada por un valor de A (lambda) sea igual a cero. 

\C-BA\ = 0 


l-0 2 (|-/?)-A(l + /?0 2 ) 
-0 2 (l-y)-A(0 2 y) 


1 

1 - A 


= 0 


Resolviendo el determinante e igualando a cero, se obtiene una ecuacion 
caracterfstica 


A 2 + 


-2 + 




i + pe 2 


A + 


1 - 


1)2 (y-S 


i + pe 2 


= o 


(17) 


En la respuesta de vibracion libre, se espera una solution de sen / cos con 
amplitud de 1; al valor de A (lambda) se puede colocar como A = Pe~ w o como 
tambien A = Pe w , el valor de P 2 = 1 para el caso de vibracion libre y se elige A = 
e w , para sustituir en la ecuacion (17). 


e 2W + 


-2 + 


1 +/?0 2 


e w + 


9' 


1 -■ 


JizM 

i + pe 2 


= o 


(18) 


La ecuacion anterior se divide para e w y se resuelve algebraicamente. 


e w + e~ w + 


9 


-2 + ■ 




1 + P9 2 


— p-iO 


— c 

1 +/J9 2 


(19) 


Consideramos la relacion de Euler que nos indica e w = cos 9 - i sin 9 y e~ w = 
cos 9 + £sin0, por consiguiente se remplaza en la ecuacion anterior obteniendo en 
funcion de senos y cosenos. 


2 cos 9 + 


-2 + 


e 2 {;+r) 

1 +P9 2 


= cos 9 


glOH) ' 

1 + /?0 2 


— i sin 9 


0 2 (>--;) 

1 + p9 2 


( 20 ) 


Evaluamos para un valor de y = 14 y obtenemos 

9 2 


2 cos 9 + 


-2 + 


1 +/39 2 


= 0 


( 21 ) 
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Ya que la amplitud de la funcion coseno va desde -1 hasta 1, se considera 
que cos 9 = -1, para determinar la siguiente inecuacion. 


e 2 


2 + 2 i + pe 2 
4 + 4/?# 2 = e 2 

9 2 (1 — 4/3) = 4 
4 

e 2 < 


1-4/? 


e < 


( 22 ) 




Se recuerda que el valor de 0 = u> At, por lo que resulta: 

4 

At < 


(23) 




Cl) 


Para finalizar se colocan los valores del metodo con y=1/2 y p=0, asi se 
determina cual es el rango que debe tener los valores que asignemos a la variable 
At para que los resultados no sean incorrectos como los especifica Aref (2016), se 
toma el criterio de Courant que es una condicion de convergencia que condiciona 
un limite del tiempo para que la simulacion no tenga errores, en este caso el Ifmite 
que tiene el tiempo esta en funcion de la frecuencia natural dando como resultado 
que: 

A 2 (24) 

At < — v ' 

Ci) 

La estabilidad del metodo de Newmark esta generada por otra inecuacion 
en funcion de la variacion del tiempo y el periodo natural de la estructura como 
mencionan Chopra A. (1995) y Clough & Penzien (1995). 

At 1 1 

< 


T n ttV2 Jy — 2/? 

Sustituyendo los valores de y=1/2 y 3=1/4 en la inecuacion anterior resulta. 

At < 0.551 T n (25) 

3.2 Procedimiento de calculo 


Para el siguiente calculo se determina como valores iniciales de 
aceleracion, velocidad y desplazamiento igual a cero, el valor de y que es 0.5, el 
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valor de (3 dependiendo del resultado de calculo puede variar entre !4, 1/6 y 1/8, que 
significa un metodo de aceleracion: constante, lineal o escalonada y At que consiste 
en una variation de tiempo con la que se desea trabajar. 

1 . Crear las matrices M, K y C del sistema. 

2. Inicializar con los valores de: 


q, q, q = 0 

3. Colocar un valor de At, y y 3, y generar constantes para facilitar el calculo 

y > 0.50; y? > 0.25(0.5 + yf (26) 

1 

- Jit 
f = -L- 

h p At 
1 

~JU 2 

1 

^~IJ 

4. Formar la matriz de rigidez efectiva 

K = K + Mf 2 + C A 

Los siguientes pasos se realizan de forma interactiva para cada piso 

5. Calculo de fuerza efectiva para el instante de tiempo i+1 . 

Fi+i = ~MJ a l+1 + M(f 0 q t + f 3 q £ ) + C(f 4 q t + f 5 At q £ ) - K q t 

6. Se determina el incremento del desplazamiento en el tiempo con la Fuerza 
efectiva y la rigidez efectiva. 


K*Aq l+1 = F i+1 
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Al resolver el sistema de ecuaciones se halla A q i+1 


7. Se determina los valores de la aceleracion, velocidad y desplazamiento 
como respuesta en el tiempo. 

Qi+i = (A 1 ) - (/o Qi) - (A <7 i) 


Qi+i = (A A(? i+1 ) + (A qd + (A Mi +1 q i ) 


Qi+i — A(?i + i + qi 


8. Los valores de aceleracion, velocidad y desplazamiento se actualizan para 
determinar el incremento en el tiempo con el calculo nuevamente. 

Ri = Qi+l 

Qi = Qi+1 
Q i = Qi+1 


4. METODO DE WILSON 

El Sistema de ecuaciones diferenciales que es usado en los problemas 
dinamicos define la ecuacion: 

Mq + Cq + Kq = M ] a(i) (27) 

Para el uso del metodo de Wilson se inicia con una valor de 0 > 1 para 
asumir una variacion lineal de la aceleracion para el paso del tiempo desde t hasta 
t+0At, pero para asegurar una estabilidad incondicional es necesario usar un valor 
de 0 > 1 .37, usualmente se emplea un valor de 0 = 1 .4. 




292 


Esteban Vasconez, Roberto Aguiar, Amjad Aref 



Para deducir la ecuacion de la recta de la aceleracion en funcion de t, se 
determina la pendiente considerando dos puntos. 


<h+9At ~ 9 i 

m — 7m : 

1 + 8 At - i 

La ecuacion de la recta por medio de dos puntos esta dada como: 

q T - i = m(j - i) 

Donde se remplaza la ecuacion (19) en (20) 

<7i+0At — <7 i . 


(28) 


Q i 


8 At 


-(r - t) 


(29) 


(30) 


Remplazando r = i + r y resolviendo de forma algebraica, con lo que para el 
intervalo de i a i + 0At se tiene: 


8i+ r ~ 8 i g ^ (Qi+BAt <7 i) 


(31) 


Integrando la ecuacion dos veces para determinar la velocidad y 
desplazamiento en funcion de r. 


8i+ x 8i "I" TCj i + 


26 At 


(Hi+dAt 8 i) 


(32) 


q t+T = q t +Tq 1+^ + ^ 


((/i+0At - <7 l ) 


(33) 
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Las ecuaciones anteriores (32) y (33) se remplaza los valores de r por 0At. 

_ • 9At r - •• > ( 34 ) 
Qi+BAt — Qi "L 2 Wi+SAt Q i) 


0 2 At 2 ,.. , (35) 

tfi+dAt ~ Qi QAtC(i H ~ vtfi+9At 2cj j) 

Se obtiene la aceleracion y velocidad de i + 0At en terminos del 
desplazamiento de i + 0At. 

La ecuacion (36) se obtiene al despejar q i+0At de la ecuacion (35) 


Qi+BAt ~ 


6 

6 2 At 2 


( Mi+BAt Qi) * 


(36) 


La ecuacion (37) se obtiene de reemplazar la ecuacion (35) en la (34) 

3 OA t„ (37) 

tfi+BAt — Q ht 2 ^ 

Para obtener la solucion del sistema y determinar los valores de aceleracion, 
velocidad y desplazamiento se usa la ecuacion de movimiento en el tiempo i + 0At 
extrapolando linealmente un vector de fuerzas, por lo que se usa lo siguiente. 


M(\i+QAt "I" CQi+BAt 4" Kqi+BAt ~ F i+9At (38) 

Para el calculo de la fuerza efectiva en funcion de i + 0At se mantiene la 
siguiente ecuacion despues de remplazar (34) y (35) en (38) y resolver 
algebraicamente: 


Fi+6At = Ft + Q(F l+ At - Ft) (39) 

Por lo que remplazando la expresion de Masa, J y aceleracion en la 
ecuacion (39) se obtiene la Fuerza 

Fi+eAt = MJ a.i + 8(MJ a i+At - M J a t ) (40) 

Ecuacion de la rigidez efectiva se obtiene despues de remplazar 
(34) y (35) en (38) y resolver algebraicamente. 


K=K+M (db) +C 


(41) 
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El calculo de desplazamiento en i + 0At, esta dada en funcion de la rigidez 
efectiva y de la fuerza efectiva, para ello se debe resolver el siguiente sistema de 
ecuaciones lineales. 


K * Qi+01 — Fi+ei ( 42 ) 

De la cual despejamos el desplazamiento. 

qi+ei = K* 1 Fi+ei ( 43 ) 

Para obtener q l+At se inicia evaluando r = At de la ecuacion (31) y se 
reemplaza la ecuacion (36) como se muestra a continuation. 

At ( 6 f ^ 6 o .. ^ 

Qi+At = Q i + 2 (Qi+ 0 At — Qi) — Qi ~ i ~ Q ij 

_ 6 f ^ 6.^3^.. (44) 

^ +At _ 0(02At 2 ) {qi+9At ~ qi) ~WAtj qi + \ 1 ~e) qt 

Para obtener q i+At y q; +At se evalua la ecuacion (34) y (35) que el valor de 
0 = 1 y se obtiene lo siguiente: 


. At 

Qi+At — Qi " I" 2 (Qi+At 4” *7 i) 


At : V 

Qi+At — Qi 4" Attjj + (Qi+At "I" 2 q i) 


4.1 Procedimiento de calculo 


(45) 

(46) 


Para iniciar con el calculo se considera el valor de 0 segun las caracteristicas 
anteriores mencionadas y corregir el metodo. 

1 . Crear las matrices M, K y C del sistema. 

2. Inicializar con los valores de: 

q, q, q = 0 

3. Colocar un valor de At y 0, y generar constantes para facilitar el calculo 

9 > 1 lineal ; 9 > 1.37 9 = 1.40 (47) 


fo = 


6 

(0 At) 2 3 * ’ 


A 


3 

61 At 5 


A = 2/1; 
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U 


fo. 
8 ' 



fs= Z j L \ fe = 1 


4. Formar la matriz de rigidez efectiva 


3 . 

8 ' 



K =K + Mf 0 + Cf 1 


Los siguientes pasos se realizan de forma interactiva para cada piso 

5. Calculo de fuerza efectiva para el instante de tiempo i+1 se calcula con la 
ecuacion (19). 


F i+ e i = MJUi+ 9(M ] a i+1 - MJ a ; ) 


Fi+ei ~ Fi+gx + M(/ 0 q + f 2 eh + 2 q i) + Cfaq + 2 4i + h q ;) 

6. Se determina el incremento del desplazamiento en el tiempo con t + 0At 
con la Fuerza efectiva y la rigidez efectiva. 

li+81 = K 1 F i+ g 1 

7. Se determina los valores de la aceleracion, velocidad y desplazamiento 
como respuesta en el tiempo. 

9i + 1 f^^Qi+dl ~ Qi) "F (Js tfi) "F if 6 Q i) 

Qi+1 = Qi+ fliAi+l + Q i) 


<Ji+ 1 = q t + At q t + fa(q i+1 + 2q t ) 


8. Los valores de aceleracion, velocidad y desplazamiento se actualizan para 
determinar el incremento en el tiempo con el calculo nuevamente. 

r?i = r?i + i 
Ri = Qi+l 
<7 i = 9i+l 


5. EJEMPLO DE APLICACION DE LOS METODOS 

Para el ejemplo de analisis se usa un portico compuesto de vigas y 
columnas, ver figuras 4 y 5, el metodo requiere una base de datos generada por un 
sismo que indica las aceleraciones en varios puntos del tiempo para este ejemplo 
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se usa el archivo rec_centro que son los resultados ante el sismo del Centro en 
1940. Las matrices de Rigidez y Amortiguamiento se obtiene usando CEINCI-LAB ; 
la matriz de masa se realiza un calculo aparte. 



Figura 4: Portico propuesto para el analisis sismico 


2 T/m 


limn 

mini 

0.3/0. 3 

0. 3/0.3 

0.3/0 .3 

0. 3/0.3 

2 T/m 

mini 

limn 

0.3/0. 3 

0. 3/0.3 

0.3/0 .3 

0. 3/0.3 

2 T/m 

m mi i 

T1TT17I 


m2^Bl.633Ts 2 /m 


ml 1 .633 T s 2 /m 


mO 1.633 Ts 2 /m 


Figura 5 Modelo de calculo del portico piano para el analisis 


El calculo de la matriz de masa se realiza usando la carga distribuida impuesta 
con respecto a la longitud del portico y la gravedad. 

w * L 


m = 


9 
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2 T /m * (4m + 4m) 
mO = ml = m2 = - 

mO = ml = m2 

Datos obtenidos del analisis en el 
caracterfsticas de la estructura como la 
matriz de masas y J son datos. 



' 3706.1 

-2062.4 

380.8 ' 


' 7.37 

-2.49 

0.01 ' 

K = 

-2062.4 

3064.3 

-1454.3 

C = 

-2.49 

6.15 

-2.45 


. 380.8 

-1454.3 

1128.1 . 


. 0.01 

-2.45 

3.52 . 


T 


1.633 

0 

0 ' 

1 

M = 

0 

1.633 

0 

.1. 


. 0 

0 

1.633. 


Para determinar los desplazamientos, velocidades y aceleraciones de la 
estructura con respecto al sismo mencionado se usa el programa 
Newmarklineal2_gama como se indica a continuation: 

[tablabn]=Newmarklineal2_gama(p,M,C,KL,J,dt,beta,gama) 

• p corresponde al nombre del archivo que contiene el acelerograma; 

• M es la matriz de masa; C es la matriz de amortiguamiento; 

• KL es la matriz de rigidez lateral del sistema; 

• J es el vector unitario para la resolution de un caso piano; 

• dt es la variacion del tiempo del acelero grama; 

• beta es el valor de la variable que depende si se considera una aceleracion 
constante 1/4, lineal 1/6 o escalonada 1/8; 

• gama es el valor de la variable de integracion en este caso es Vz\ 

• tablabn son los valores resultantes de la interaction que entrega el programa 
en una tabla, consta del numero de pisos, aceleraciones maxi mas, 
velocidades maximas y desplazamientos maximos de cada uno de estos. 

Para estabilizar el metodo de Newmark se debe considerar que la variacion 
del tiempo At debe ser menor a 2/w como se especifica en calculos anteriores con 
el criterio de Courant, la frecuencia natural de la estructura se determinan con 
CEINCI-LAB para cada piso. 


9.8 m/s 2 
= 1.633 Ts 2 /m 

programa CEINCI-LAB, nos indica las 
matriz de rigidez y amortiguamiento, la 


(0 = 


10.5534' 

34.1726 

.59.6545. 


Por lo tanto el valor de 2/oj para cada piso es: 
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0.1895' 


‘At < 0.1895' 

0.0585 


At < 0.0585 

-0.0335. 


.At < 0.0335. 


Para comparar se usa la ecuacion de Chopra (1995) y Clough & Penzien 
(1995), conociendo el periodo natural de la estructura calculada con CEINCI-LAB. 


T = 


0.5954' 

0.1839 

.0.1053. 


Por lo tanto los valores de 0.551 * T periodo natural son: 



0.328 


At 

< 0.328' 

0.5517 = 

0.101 


At 

< 0.101 


.0.058. 


.At 

< 0.058. 


Considerando los criterios anteriores en funcion de la frecuencia natural y 
el periodo natural de la estructura para estabilizar el metodo de Newmark, y 
conociendo que el archivo del acelerograma tiene un dt = 0.01, este valor de la 
variacion del tiempo es el optimo para no generar errores en el calculo, se considera 
P=1/4 y y=1/2 para una aceleracion constante. Se obtienen los siguientes resultados 
y graficas. 

Tabla 1 Aceleraciones maximas, desplazamientos maximos y velocidades 


maxi mas de cada piso con p=1/4 y y=1/2 


Nro. Piso 

Aceleracion 

[m/s 2 ] 

Desplazamien to 

[m] 

Velocidad 

[m/s] 

1 

5.506 

0.0237 

0.2348 

2 

9.374 

0.0540 

0.5441 

3 

10.602 

0.0734 

0.7970 
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Figura 6: Desplazamiento en el tiempo del piso 3 con (3=1/4 y y=1/2 


En el siguiente caso al metodo de analisis de Newmark se considera los 
valores de y = 1 /? > 0.25(0.5 + y) 2 ; entonces (3= 0.5625 para una aceleracion lineal 
en el mismo programa dando como resultado la siguiente tabla y grafico. 


Tabla 2 Aceleraciones maximas, desplazamientos maximos y velocidades 


maximas de cada piso con |3=0.5624 y y=1 


Nro. Piso 

Aceleracion 

[m/s 2 ] 

Desplazamiento 

[m] 

Velocidad 

[m/s] 

1 

4.8308 

0.0216 

0.1966 

2 

8.6545 

0.0493 

0.4755 

3 

10.1511 

0.0672 

0.7000 
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Figura 7 Desplazamiento en el tiempo del piso 3 con 3= 0.5625 y y=1 


Con el siguiente grafico se pretende comparar los resultados en el 
desplazamiento del metodo de Newmark en un rango considerable del tiempo para 
observar la diferencia al momento de cambiar los valores de beta y gama en el 
sistema de ecuaciones, el grafico muestra que al aumentar el valor de gama y por 
lo tanto el valor de beta se genera una disminucion de desplazamiento en los 
resultados a traves del tiempo. 



Figura 8 Desplazamiento en el tiempo M. Newmark con diferentes valores de 

y yp- 
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De otra forma para calcular las aceleraciones, velocidades y desplazamientos 
con el metodo de Wilson es necesario el uso del programa metodo_wilson que se 
indica a continuation. 


[tablamw]=metodo_wilson(p,M,C,KL,dt,teta,J) 

• p corresponde al nombre del archivo que contiene el acelerograma; 

• M es la matriz de masa; 

• C es la matriz de amortiguamiento; 

• KL es la matriz de rigidez lateral del sistema; 

• J es el vector unitario para la resolution de un caso piano; 

• dt es la variation del tiempo del acelero grama; teta es el valor de la variable 
que ayuda a estabilizar el metodo de calculo, debe ser mayor a 1 .37 pero se 
considera el valor de 1 .4; 

• tablamw son los valores resultantes de la interaction que entrega el 
programa en una tabla, consta del numero de pisos, aceleraciones maximas, 
velocidades maximas y desplazamientos maximos de cada uno de estos. 


El archivo del acelerograma tiene un dt = 0.01, se considera 0=1.4. Se 
obtienen los siguientes resultados y graficas. 

Tabla 3: Aceleraciones maximas, desplazamientos maximos y velocidades 


maximas de cada | 

piso usando 0=1 .4. 

Nro. Piso 

Aceleracion 

[m/s 2 ] 

Desplazamien to 

[m] 

Velocidad 

[m/s] 

1 

5.5183 

0.0238 

0.2287 

2 

9.3054 

0.054 

0.5381 

3 

10.5837 

0.0734 

0.7956 
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Figura 9 Desplazamiento en el tiempo del piso 3 con 0=1 .4. 

6. RESULTADOS 

La solution con cada metodo usando el sismo del Centro han sido similares 
en este caso generando resultados bastante parecidos tanto en las aceleraciones, 
las velocidades y los desplazamientos como se muestra en la tabla 3 


Tabla 3 Resumen de resultados. 



Newmark 

Newmark 

Wilson 


P=0.25 y=0.5 

P=0.5625 y=l 

0 = 1/4 

Nro. Piso 

Aceleracion [m/s 2 ] 

Aceleracion [m/s 2 ] 

Aceleracion [m/s 2 ] 

1 

5.5060 

4.8308 

5.5183 

2 

9.3740 

8.6545 

9.3054 

3 

10.6020 

10.1511 

10.5837 


Desplazamiento [m] 

Desplazamiento [m] 

Desplazamiento [m] 

1 

0.0237 

0.0216 

0.0238 

2 

0.0540 

0.0493 

0.0540 

3 

0.0734 

0.0672 

0.0734 


Velocidad [m/s] 

Velocidad [m/s] 

Velocidad [m/s] 

1 

0.2348 

0.1966 

0.2287 

2 

0.5441 

0.4755 

0.5381 

3 

0.7970 

0.7000 

0.7956 
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Los dos metodos con las consideraciones previas nos arrojan resultados 
bastante similares en el desplazamiento. 




Figura 10: Esquema del desplazamiento en el portico de analisis. 


En el siguiente grafico tornado en un rango de 25 segundos se puede 
observar la diferencia entre los metodos y cambiando las variables gama y beta. 



Figura 11 Desplazamiento en el tiempo en 8 s 


6. CONCLUSIONES 


• Obtener resultados del comportamiento de la estructura con respecto a 
una solicitacion dinamica implica conocer que tipo de dano puede sufrir 
dependiendo de su material. 

• Los metodos fueron deducidos por integration y se generan por medio 
de interacciones, son muy sensibles a las variables iniciales como (3, y, 
0 y At. 
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• Se determino el valor limite para la variacion del tiempo At que depende 
del criterio de Courant, el valor no puede ser mayor a 2/frecuencia o 
tampoco a 0.551 *Periodo como expresan dos autores para evitar 
errores en los resultados en el metodo de Newmark. 

• Dado la aceleracion de piso inicial con el Sismo del Centro se 
determinaron valores similares en los dos metodos que van ser usados 
para conocer el comportamiento de la estructura. 

• Los resultados al colocar un gama y un beta mas elevados en el sistema 
de ecuaciones usado en el metodo beta de Newmark son menores tanto 
en aceleracion, velocidad y desplazamiento en cada piso por lo que 
concluye el uso de estas variables si existe amortiguamiento en la 
estructura. 
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